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Abstract. The main objective of this article is to study the effect of spherical geometry on 
dynamic transitions and pattern formation for the Rayleigh-Benard convection. The study is 
mainly motivated by the importance of spherical geometry and convection in geophysical flows. It 
is shown in particular that the system always undergoes a continuous (Type-I) transition to a 2lc- 
dimensional sphere S^'", where Ic is the critical wave length corresponding to the critical Rayleigh 
number. Furthermore, it has shown in [12] that it is critical to add nonisotropic turbulent friction 
terms in the momentum equation to capture the large-scale atmospheric and oceanic circulation 
patterns. We show in particular that the system with turbulent friction terms added undergoes 
the same type of dynamic transition, and obtain an explicit formula linking the critical wave 
number (pattern selection), the aspect ratio, and the ratio between the horizontal and vertical 
turbulent friction coefficients. 
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1. Introduction 

The main objective of this article is to study the Benard convection of fluids on a spherical shells. 
The main motivations are two-fold. First, convection plays an important role in the large scale 
atmospheric and oceanic circulations. The phenomena of the atmosphere and ocean are extremely 
rich in their organization and complexity, and involve a broad range of temporal and spatial scales 
[3] . To avoid to deal with the atmosphere and ocean in all their complexity, the main focus of this 
article is to study the Benard convection on spherical shells to capture the main effects caused by 
the spherical geometry. Second, over the years, the Rayleigh-Benard convection problem, together 
with the Taylor problem, has become one of the paradigms for studying nonequilibrium phase 
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transitions and pattern formation in nonlinear sciences; see among others Busse [1] , Chandrasekhar 
[2], Cross [4], Ma and Wang [13], and the references therein. 

For the Euchdean domain case, a systematic study on dynamic transition of the Rayleigh-Benard 
convection is carried out in Ma and Wang [9, 11]. The main result in this direction is that the system 
always undergoes a Type-I (continuous) transition as the instability driving mechanism, namely 
Rayleigh number, crosses a critical threshold i?c, thanks to the symmetry of the linear operator, 
properties of the nonlinearity and asymptotic stability of the basic state at the critical threshold. 
Moreover, the system has a bifurcated attractor which is an {m. — 1) -dimensional homological sphere 
where m is the number of critical eigenvalues of the linear operator. 

In this article, we consider a fluid which occupies a layer of spherical shell modeled by the 
Boussinesq equations. Namely, the density function is treated as a constant except in the equation 
of state and in the buoyancy term. For simplicity, the spatial geometry of the domain is taken as 
the product of a two-dimensional sphere with radius a and an interval: x (0,/i), where h is 
the height of the fluid layer. This is mainly motivated by the fact that aspect ratio for the large 
scale atmosphere and ocean is small; see among others [7, 8, 14]. 

As in the Euclidean case, in the spherical shell case, we show that the system is always undergoes 
a Typc-I (continuous) dynamic transition to a 2Zc-dimensional homological sphere Sjj, where Ic is 
the critical wave number determined by ^ 

(1.1) 

An important part of this article is to show that 2?c-dimensional homological sphere Eij is in fact 
homeomorphic to a 2Zc-dimensional sphere S'^^". This is achieved by carrying out the reduction of 
the original Boussinesq equations to the center manifold generated by the first 2lc + 1 modes. The 
main difficulty comes then from the nontrivial nonlinear interactions of the spherical harmonics- 
both scalar and vectorial. 

The remaining part of this article is devoted to the introduction of the turbulent friction terms 
in the Boussinesq equations. This is crucial for capturing the large scale convection patterns in the 
large scale geophysical flows, as demonstrated in thermohaline circulation studies by Ma and Wang 
[12]. The friction terms are added based on the convection scale considerations as introduced in 
[12], and are given by {aou,aiw); see (6.1). Under the assumption that 

1 <C (70 <C CTl , 

the corresponding critical wave number Ic is approximately given by 
(1.2) 

As we know, for the large-scale atmospheric and oceanic circulations, the critical wave number 
is often relatively small, e.g. Ic = Q for the global Walker circulation over the tropics. Hence the 
small aspect ratio h/a is then balanced by the smallness of the ratio ao/ai between the horizontal 
and vertical turbulent friction coefficients, leading to moderated size of the critical wave numbers. 

It is worth mentioning that with the framework given in this article, one can study the humid- 
ity/salinity effects of the large scale atmospheric and oceanic flows, as well as the rotational effect 
(Coriolis force) of the earth. We shall report this elsewhere. 



lc{lc + l) = 



h 
a 



Icilc + 1) ~ Y 



CTl 



The integer is to make the left-hand side the closet to the right-hand side. 
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The article is organized as follows. Section 2 introduces the Boussinesq equations on the spherical 
shell, and the linear stability analysis is carried in Section 3. Sections 3 and 4 are on dynamic 
transitions and the structure of the bifurcated attractor. Section 6 demonstrates the need to 
introduction of convection scales, and studies the dynamics transitions for the Boussinesq equations 
with the added friction terms. 

2. Boussinesq Equations 

We start with the prototype of problem, the Benard convection problem for a layer of fluid in a 
spherical shell. We consider the case where the fluid can be modeled by the Boussinesq equations. 
Namely, the density is considered as constant except in the buoyancy term and in the equation of 
state. The equations are given as follows: 

ut + {u- V3)u - z/Aau + Po^VsP = -gt{l - a{T - %)), 
(2.1) Tt + {u- V3)T - kAsT = 0, 

divsu = 0. 

Here the unknown functions are velocity fleld u = {ui, U2, Us), the temperature function T, and the 
pressure function p. Other parameters include: ly is the kinematic viscosity, k = (0, 0, 1) is the 
unit vector in the vertical direction, g is the gravity constant, a is coefficient of thermal expansion 
of the fluid, Tq is temperature on the lower surface, po is the density,«; is he thermal diffusivity.For 
simplicity, we consider the spherical domain as 

nd = slx io,h), 

where is the two dimensional sphere with radius a, and h is the height of the fluid layer. The 
motion for this type of geometry is due to the study of large scale geophysical fluid flows, where 
the aspect ratio, and the thin layer of spherical shell can be approximated by fid given here. 
It's easy to get the basic solution: 

{u,T) = iO,%-pX3), 

p = po-pog{x3 + -^^l)- 

where 

and Ti is temperature on the upper surface. 

Mathematically speaking, it is convenient to introduce the nondimensional form of the equation. 
For this purpose, let 

{x,t) = {hx', — ), 

K 

{u,T) = {^,Ph^ + Jh-l3hx',), 

P = Po-f^+Po - gpoihxs + (X3) ), 
In addition, the Rayleigh number R and the Prandtl number Pr are deflned by 
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The nondimensional domain is 



n = s^x (0, 1), 



and it's reasonable to use the spherical coordinates. Let {u,w) be the 3 dimensional velocity, u be 
the horizontal velocity field, and w be the vertical velocity. Also the following notations div, V are 
differential operators in the horizontal direction. Then the system (2.1) can be written as: 



(2.2) 



+ V„w + u,^ + Vp) - (A + -^)u = 0, 



Tt + V„T + w 



dw 



divu ■ 



dw 



dz 
0, 



^^«'-(A + ^)^ = 0, 



supplemented with the following free-slip boundary condition: 



(2.3) 



du 

w = 0, T = 0, — = at z = 0,1. 
oz 



It is worth mentioning that we can also study other physically sound boundary conditions such as 
the non-slip boundary or the periodic boundary conditions. 

The following are some differential operators at spherical coordinates on S^: 



„^ Idf 1 df 

r oO r sm dip 



curlT 



1 df^ 

r sin 9 dcp 



ldf_ 



Uu, dv0 



dv0 

86 ' smO dip 
V„T= - (we — 



sm.9 dip 



AT: 



r 

1 1 
sin 9 



dT 



sin 6 dip ' 
d , . dT, 1 d'^T 



2 cos 9 duu, 
Au = {Au0 ^ 



U0 



2 oln^ Q(_p 



sm 



sin^ 



-)ee + (Au^ + 



2 cos 9 du0 
r'^sm^9 dip 



sin^ ( 



where u, v are vectors, and T is scale. 

For the problem (2.2) with (2.3), we define the function spaces as 



(2.4) 



H = {{u,w,T) e L^(0)^| divu = 0;{w,T) =0, at ^ = 0, 1}, 

Hi = {{u,w,T) e H^{Q)^\ divu = 0;{w,T) =0,^ =Oat ^ = 0,1}. 

dz 



Let Lx'- H, and G : Hi ^ H he defined by 
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Lx6 = P 



(2.5) 



02 

Pr(A+ ^)«; + PrAT 

' du ' 

VuU + W — 

oz 

„ dw 

VuW + W — 

oz 

„ dT 
VuT + w— 
oz ■ 



for any cf) = {u, w, T) <E H^. Hero A = ^/R, and P : L^(ri)^ H is the Lcray projection. 
Then the equations (2.2) with (2.3) can be rewritten in the following operator form: 

(2.6) 



^ = L,<^ + G(<^), 



for t = 0. 



Here (l)o G H is the initial condition. It is then classical to prove the existence and properties of 
solutions; see [6, 8]; we omit the details here. 



3. Linear Analysis and Principle of Exchange of Stability 

To study the dynamic transitions of the Benad convection problem (2.2) on the spherical shells, 
we need to consider the linear eigenvalue problem: 

(3.1) = ^(A)</>, 

which is equivalent to 



Pr(A + -Vp = /3(A)u, 



(3.2) Pr(A+^^). + PrAT-|=/3(A)«,, 

{A + -^)T + Xw = l3{X)T, 
div3(u,w) = 0, 

supplemented with the boundary condition (2.3). 

When {u, w, T) doesn't depend on vertical z, then we have (w, T) = 0, and u = cmlYi^m, so the 
eigenvalue and eigenfunctions are as follows: 

PiQ = — Pra?, 

(3.3) ^° , V 

*(r„o= (curir,„,0,0). 
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Thanks to the spherical geometry of the domain, the above eigenvalue problem can solved using 
the following separation of variables: 

u = Vf{0,ip)H'iz), 

(3.4) w = a\fi0,ip)H{z), 

T = f{9,^)e{z). 

By (3.4), we infer from (3.2) and (6.2) that 

(3.5) - A/ = aV, 

(3.6) Pii-a'^VfH' + VfH'") -Vp = pVfH', 

(3.7) Pv{-a^fH + a^fH" + A/G) - |^ = pa^fH, 

(3.8) - + /e" + Xa^fH = 0fQ, 

(3.9) {e, H, H") =0 at z = 0,1. 

For the eigenvalue problem (3.5) on the 2D sphere, the only possible candidates for the eigenfunction 

and eigenvalues are given by 



(3.10) f^Yi^i9,ip) 

(3.11) af = 



dz 

PrH"" - (2Pra2 + j3)H" + (a^Pr + l3a'^)H 



To eliminate the pressure function p, we do M - V(3.7): 

(3.12) 6= 

Then we can consider two cases. First, for the case where Z = 0, it is easy to see that the 
eigenvalues and eigenvectors of Lx are: 

(3.13) MOn 

*oon = (0,0,smn7r^). 
When I 7^ 0, plug (3.12) into equation (3.8), then H, Q are as follows: 

(3.14) Hn = sinriTT^;, 

(3.15) Qn = binsmmrz. 

Then plugging (3.14) and (3.15) into (3.12) and (3.8), we obtain that 
where 

2 2 2 I 2 2 2 I "I" 1) 

Prom (3.16), the eigenvalue /? solves the following quadratic equation: 

+ (1 + Pr)7?„/3 + [it - ^)Pr = 

Tin 



(3-16) hn = a , 2 
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we get two types of eigenvalues and eigenvectors: 



^ lln 



(3.17) 



= ("'^cos 71772; Vyjm, afYimSinmrz, —r — '-^YimSinmrz), 

Pin = JhlfJl + Pr) - (7f„(l + Pr)^ - 4Pr(7f„ - ^))^/^], 
^ rin 



2 



*/mn = ('^■^COS 71772; VYim, afl^m sin 71772, ' „ sin 71772) . 

So for problem (2.2) with boundary condition (2.3), the eigenfunctions for the corresponding 
linear operator as in (3.2) are consist of three group, which is shown in (3.3), (3.13), and (3.17). 
Now we define the critical Rayleigh number Rc = by 

(3.18) Ac = min^^ — = min ^^—^ — . 

i.n ai I ai 

It is easy to sec that the minimum in the above formula is Ac is achieved at of . By definition, we 
have (see (1.2)): 

'h] 



(3.19) 



2 2 

lc{lc + l) = Y- 



As I can only take integer values, it is possible for the minimum in (3.18) is achieved at two 
consecutive Z's. For simplicity and without loss of generality, in this article, we only consider the 
case where the minimum is achieved at ai^ with an integer Ic- 

X 3^/3772 /c(?e + l) 

(3.20) Ac = = — - — with ai^ = 5 , 

and the critical Rayleigh number is 

2777* 

(3.21) R^ = xl = -— = 657.5. 

This value is consistent with the value developed classically; see among many others [2, 5]. 
By definition, the following principle of exchange of stability (PES) holds true: 

r <0 if A < Ac, 

PtiWl =0 if A < Ac, 
[ >0 if A > Ac, 

/3;+(Ac)<0 V(/,7i)^(Zc,l), Z,7i>0 

;3-(Ac)<0 y {l,n), l,n>0, 

/3io(Ac) < V / > 0, 

/3on(Ac) < V 71 > 0. 
The critical eigenspace is then given by 

Ic 

(3.22) ^1 = { Xm^tmi I ^-m = (-l)™^;;,, Xm & C, \m\ < Q 
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with dimension dimi^i =2lc + 1- 

4. Dynamical Transitions 

As mentioned in the introduction, the hnear operator Lx is symmetric, leading to the real 
eigenvalues (3f^^^, /3o,n- With the PES at our disposal, we can obtain immediately that the Benard 
convection system will undergoes a dynamic transition, and the following theorem shows that the 

transition is always Typc-1. 

Theorem 4.1. The Benard convection problem (2.6) always undergoes a Type-I dynamic transition 
as the Rayleigh number R crosses the critical value Re- In particular, the following assertions hold 
true: 

(1) When R < Rc, the basic motionless state with linear temperature profile cj) = {u,w,T) = 
is globally asymptotically stable. 

(2) As R > Rc, but near Rc, the problem bifurcates to an attractor which is an 2lc- 
dim,ensional homological sphere, where the integer Ic is defined by (3.20). 

(3) For any (p = {u, w, T) G S^j, 

Ic Ic 

X] ^™*tml+0( ^ \Xm\), X-m = {-'^)'^X*^, Xm & C, \m\ < Ic 

— _1 ™ — ; 



(4) The attractor Yir attracts H\T, where T is the stable manifold of tj) = with codimension 
2lc + l. 

(5) The attractor Ei{ contains at least one steady state solution of the problem. 

Proof. As in [9] , we can show that the basic motionless state is (globally) asymptotically stable the 
critical eigenvalue A = Ae- By the attractor bifurcation theorem in [10] and the linear analysis in 
the previous subsection, we obtain immediately Assertions (l)-(4) in the theorem. 

It suffices then to show the existence of at least one steady state in the bifurcated attractor Sjj. 
To this end, the steady state problem for (2.6) can be written as 

(4.1) 

where A - 



Bx 



-A<j)-Bx^^Gi<j)), 
Lx and the linear operator A and Bx are defined by 



(4.2) 



A(I) = P 



52 

Pr(A+^)u 



Pr(A- 



Bx<P = P 





PrAT 
Xw 



\/(j>={u,w,T) e Hi. 



Obviously, the operator A is invertible, and (4.1) can be written as 
(4.3) [id + {-A)-'Bx]cl^ = i^A)-'G{q^). 

Then it is classical to show that the linear operator id + {—A)~^Bx ■ H H is a completely 
continuous field, and {—A)~^G : H ^ H is a. compact operator. Hence by the Krasnoselski 
bifurcation theorem (see e.g. Theorem 1.10 in [10]), (4.3) has a nontrivial steady state bifurcation 
at Rc, as the critical eigenvalue Rc has odd multiplicity 2lc + 1, which is the dimension of the 
eigenspace Ei . The proof is complete. □ 
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5. Structure and Patterns of the Benard Convection 

The above theorem shows that the system always undergoes a dynamic transition as R crosses 
the critical Rayleigh number Rc, and the transition states occupy a set T,fj, homological to 5^'"". 

The main objective of this section is to derive more detailed structure of this bifurcated attractor. 
In particular, we shall give a strategy to show that S/j is in fact homeomorphic to 5^'"=. 

The approach for this purpose is reduce the governing partial differential equations to the center 
manifold generated by the unstable modes for R> R^ and near R^. As the general case is tedious 
and time consuming, we present here only two special cases where critical wave number = 1 or 
2. The other cases can be studied in the same fashion. 

Theorem 5.1. Consider the Benard convection in the case where /c = 1, which is equivalent in the 
dimensional form to 

h TT 

a ^ 2' 

where a is the horizontal length scale and h is the vertical scale. Then the bifurcated attractor Ttn is 
homeomorphic to S'^, consisting of only degenerate steady states. Furthermore, is approximated 
by 

1 

m— — 1 

where qi is given by 

37r3(17787 + 355912Pr- 669713Pr^ + 387787Pr^) 



qi = 

Namely, for any cf) G S^, 



400000Pr(353 - 625Pr+ 353P7^) 



1 



m=— 1 

a;^ + |a:_i|2 + |xi|2 =/3+/gi, 
x-m = {-irx*^, \m\<l. 
Theorem 5.2. For the case where lc = 2, equivalently 

h TT 

a ^ 2^3' 

the bifurcated attractor J^r is homeomorphic to S^, which consists of at least an subset of 
degenerate steady states. Furthermore, when Pr = 1, is approximated by 

2 2 

S^~{E*L2 I = (-!)"<, E l^-l' = ^2+1/92}, 

m=— 2 m=— 2 

where q2 is given by 

_ 22914057r3 
~ 214754176' 

The proof of these two theorems is based on the reduction of the original Boussinesq equations 

to the center manifold generated by the first unstable modes, and to study the dynamics of the 
reduced system. As mentioned, the same method applies to other cases as well. 
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Proof of Theorem 5.1. We proceed in several steps as follows. 

Step 1. General Center Manifold Reduction Stratege. More precisely, let 

Ei = \x= J2 2^™*itmi I a;-m = (-l)"a;;„ x„, e C, |m| < le\ , 

I m=-lc J 

The main idea of the center manifold function is to seek a function which maps a neighborhood 
of the original in Ei into E2 given by 

(5.1) i>(a;,A)= ^ Vimo'^imO + Y, 

1>1, \m\<l n>l 

(i,n)5!t(<,4,i,n>l |m|<i i,n>l|m|<i 

By the approximation formula derived in [13], the center manifold function $ can be approximated 

by 

(5.2) -Lx^x, A) = P2G{x) + 0(2), 
where P2 : -ff — ^ -^2 is the canonical projection, and 

0(2) = o(|Mp) + 0(1/3+1 (A)|.|N|^ 

The first term in the right-hand side of (5.2) involves products of the eigcnfunctions ^'+„ji, which 
are defined in terms of the spherical harmonics Yim with I < Ic. Hence by the properties of the 
3-j symbols for the products of spherical harmonics, we only have to calculate the following center 
manifold function coefficients and all the others are zero or of higher-order: 

(1) For 1 < Z < 2lc, \m\ < I, 

(5.3) yimo = ^-M-jT 112— + 0(2); 

(2) For n = 2, 

/t, ,N < G(a:),^oo2 > ^ 
(5-4) yoo2 = — 77^ + 0(2); 



„5o2||* 



002 



(3) For ; ^ 1 < ^ < 2Zc, \m\ < and n = or 2, 

r^ln 1 1 imn 1 1 

(4) For 1 < Z < 2/c, \m\ < Z, and n = or 2, 



(5-6) yr^„ = - :-7;.-'T. 



< G(X)^;„„ > 

An 1 1 ^(mn I 



Wc note that hand-calculation of the right had sides of above formulas arc still difficult if not impos- 
sible. Fortunately, one can write a simple mathematica program to perform symbolic calculations 
to derive precise formulas these coefficients. 

Finally, the reduced equations are given as follows: 

dx^ + , <G(X + $(.T,A)),^+„, > 
(5.7) -^=PtiWxm + — + -+ ^ for |m|<Ze. 
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Again, to carry out the calculation for the nonlinear interaction in the right-hand side of this 
reduced equation, we need to use mathematica as well. 

Step 2. Approximation of the center manifold function for /c = 1- In this case, the 
first eigenvalue of (3.1) with (3.4) is 

PtiW = ^[-97r^(l + Pr) + V8l7r4(-l + Pr)2+48PrA2]. 
Using Mathematica, we infer from (5.3)-(5.6) that: 

P22 P22 

vt-ui^^ = — ^\/2rffa;-ia;o, ^^-12(2:, A) = —^V2d^x-iXo, 

P22 P22 



(5.8) 



where 



2/202(2;, A) = —^^dt{xl+x-ixi), ^202(2;, A) = 3^y|c^i (a^o +a;-ia;i), 

2/212(2;, A) = —^V2d'^xoXi, y^^^{x,X) = —^V2d^xoXi, 

^,"22 ~P22 

1 2 1 2 

2/222(2;, A) = — 2-11 2/222(2;, A) = — 2;i, 

^p22 ^P22 

2/002 = -^VSir'^ixl - 2x-iXi), 

2/;t„ = 0, when(Z,n)9^(2,2), />0. 

"3" 7r^/2(121- 121Pr + A)(187- 121Pr + A) 



dt = S\ 



10 1936[-1493 - 1331Pr^ - 11^ + Pr(2500 + llA)] ' 



o A3", 5/2. (-121 + 121Pr + A)(-187+121Pr + A) 



10 ' ' 1936[1493 + 1331Pr^ - llA + Pr(-2500 + llA)] 
/32t2 = ^^'(-121-121Pr + A), 

/32";2 = -^'^'(121 + 121Pr + ^), 



A = VTlV 1331 - 2338Pr + 1331Pr2. 

Hence wc obtain immediately the following second-order approximation of the center manifold 
function for A > Ag: 

(5-9) $(X, A) =2/2t22*t22 + 2/2^12*^-12 + 2/2^2*202 

+ 2/212*212 + 2/222*222 

+ 2/2-22*2-22 + 2/2-12*2-12 + 2/202*202 
+ 2/212*212 + 2/222*222 + 2/002*002 + 0(2), 

with the coefficients given by (5.8). 
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Step 3. Reduced equations. With the approximation of the center manifold function above, 
we deduce immediately following reduced equations: 

dx-i 



dt 



- qix-i{xl - 2x-ixi) + o(3) 



(5.10) ^ = /3+ xo - qiXo{xl - 2x-ixi) + o(3), 

^ = ^^^xi - qixi{xl - 2x-ixi) + o(3). 



where 



dt 

37r3 (17787 + 355912Pr - 669713Pr2 + 387787Pr3) 
400000Pr(353 - 625Pr + 353Pr^) 
o(3)=o(||x||3)+0(|/3+J.||x||3). 

As the eigenfunction are complex, we can set 
By definition (see (3.22)), we have 

= Zm = {-ir+'^Z_m, m > 0. 

Hence we obtain the following reduced equations: 

^ =0tiyi - Hyii4 + yl + zl) + 0(3), 
(5.11) ^ =13+ xo - q,xo{x'o + 2/? + ) + o{3), 

=/3^iZi - qizi{xl +yl + zl) + o(3). 

Step 4. Let g be the cubic terms in the right hand side of the reduced equation. Then, we have 

-{9,{yuXo,zi)) < qi{xl + yl + zif . 

Hence by the attractor bifurcation theorem in [13], the bifurcated attractor Y^r is homeomorphic 
to S^. 

Also, by Theorem 4.1, there is at least one steady state solution in Eij. In addition, due to the 

spherical geometry, the original Boussincsq equations have an S'^-symmctry. Hence, this steady 
state solution generates an ^^-set of steady states. Therefore, S/j consists precisely steady state 
solutions. 

Finally, by setting the right hand-side of the above reduced equations (5.11) equal to zero, we 
obtain that S/j is approximated by 

{(2/i,a;o,zi) gM^ | xl + yl + zl = ^ti/qi}- 
The proof is complete. □ 

Proof of Theorem 5.2. We note that the unstable space Ex is now given by 

Ei = \x= ^ a;„*2mi I x-m = {-l)'^xl„ Xm & C, |m| <lc\, 

I m=-2 J 
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and we still use for brevity x to denote 



2 

m=-2 



As in the previous proof, the approximation of the center manifold function is then given 



(5.12) <^{X, A) =y}_22'^t-22 + yt-12'^t-12 + y}o2'^to2 + 1/^12*212 + 2/2^22*^22 

+ 2/4142*4-42 + 2/4'-32*4-32 + 2/4^-2,2*4-22 + 2/4"- 12*4- 12 + 2/402*402 
+ 2/412*412 + 2/422*422 + 2/4*^3,2*432 + 2^442*442 

+ 2/2-22*2-22 + 2/2-12*2-12 + 2/2,0,2*202 + 2/212*212 + 2/222*222 

+ 2/4-42*4-42 + 2/4-32*4-32 + 2/4-22*4-22 + 2/4-12*4-12 + 2/402*402 

+ 2/^12*412 + 2/422*422 + i'432*432 + 2/442*442 + ^002*002 + o(2), 



14 WANG AND YANG 

where the non vanishing coefBcients in the above center manifold function are as follows: 



2/2^22(^1 ^c) — 
y2-12ix,Xc) = 
2/202 (•^' -^c) = 

^212 (a;, Ac) = 

2/222 C-^' Ac) = 

2/4^42 (•^' Ac) = 
2/4-32(2;, Ac) = 



^2 ' 
2Ci {X-iXq — ^f&X-2Xi) 

^2 ' 

2cI{xq - - 2x^2X2) 

2cf (xqXi - \/QX-iX2) 
^2 ' 

Ci{^/Qx\ — 4x0X2) 

^2 



"2 „2 



s/2i 



.+ 



— X-2X-I, 



2/4-22(2;, Ac) = -^{V2x\ + V3X-2Xo), 

~Pi2 



Pi2 
+ 



yt-i2ix, Ac) 



-Pt2 



{VQx-iXq + a;_2a;i), 



2/402(2;, Ac) = -^^{"ixl+Ax-iXi +X-2X2), 

~Pi2 

2/412(2;, Ac) = — -:p-(V6xoXi + X-1X2), 

~PA2 

2/422(2;, Ac) = -—^{V2xl + ^8x0x2), 

~P42 

+ V2c+ 
2/432(2;, Ac) = — r:pa;ia;2, 

~P42 



2/442(2;, Ac) 



"2 ^2 



_|_ 2:2, 



2/2-22(2;, Ac) - 
2/^-12(2;, Ac) = 
2/202 (2;, Ac) = 
2/^12(2;, Ac) = 
2/222(2;, Ac) = 
2/4-42(2;, Ac) = 
2/4-32(2;, Ac) = 
2/4-22(2;, Ac) = 
2/4-12(2;, Ac) = 
2/102(2;, Ac) = 
2/412(2;, Ac) = 
2/422(2;, Ac) = 
2/^32(2;, Ac) = 
2/442(2;, Ac) = 



Ci (v^x?_i - 4x_2a;o 



^22 



2Ci {X-iXq — \/6X-2Xi) 



'22 



{x\ - X-iXi - 2x^2X2) 

2a[{xoXi - VQX-1X2) 
^^22 ' 
Ci{\/6xl — AXQX2) 



22 



2 J2 



-Pa 
\f2c 



-P 



-X-2X-X, 



42 



-P 



-(V^X^l + \fZX-2XQ), 

-{■\f^X-\XQ +x_2a;i), 



42 



42 



-A 

~^2 
V2C2 



-{2>x1 + Ax-iXi + a;-2a;2). 



{VQxqXi + X-1X2), 

{V2xl + V^XoX2), 



X1X2, 



42 



2 ^2 



'42 



-P 



X21 



'42 



2/002 = -^y37r*(a;o - 2a;_ia;i + 2a;_2a;2), 
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where 



Co — 



(5.13) 



Co = 



y57r^/^(9 - 9Pr + ■B)(15 - 9Pr + B) 
'672[-29 - 27Pr^ - 3B + Pr(52 + 3B)] ' 

'i^f±Tr5/2(_i 153 \ 

'^V 14" \ ^ (289-289Pr+^/l7C) ^ 
'■ IQfifl I 20655 N ' 

^ -r (289-289Pr+^/l7C)2 ^ 

377^/2(442 - 289Pr + \/l7C)(289 - 289Pr + VvfC) 
16184(-11041 - 9826Pr2 - 34VV7C + Pr(18437+ 34\/l7C)) ' 

V5tt^/^{-9 + 9Pr + B)(-15 + 9Pr + B) 
672(29 + 27Pr2 -W + Pr(-52 + 3B)) ' 

3,/X7r5/2/ 1 _| 153 X 

14" V ^ ^ (-289+289Pr+vT7C) ^ 

Loyjy i-r (_289+289Pr+\/T7C)2 

-37r5/2(_442 + 289Pr + \/T7C)(-289 + 289Pr + Vl7C) 



Here 



16184(11041 + 9826Pr'' - 34V17C + Pr(-18437 + 34vT7C)) 



B = VSI- 150Pr + 81Pr2 



C = V4913^^8611Pr + 4913P?, 

and f32. 2, 1^4.2' 1^2.2' 1^4.2 are given by (3.17). 

For simplicity, we only consider the case where Pr = 1. Then the reduced equation are as 
following: 



dt 
dx-i 

dt 



=^21^-2 - q2X-2{xl - 2x-ixi + 2X-2X2) + 0(3), 
=j3}^x_i - q2X_i{xl - 2x_iXi + 2x^2X2) + o(3), 



(5-14) ^ =PtiXo - q2Xo{xl - 2X-XXX + 2x^2X2) + o(3), 



dt 
dx\ 

~dr 

dX2 

~dt 



where 



Let 



Then 



=1^21X1 - q2Xi{xl - 2x-iXi + 2X-2X2) + 0(3), 
=l32iX2 - q2X2{xl - 2x-iXi + 2X-2X2) + 0(3), 



_ 22914057r3 
~~ 214754176 ■ 



Vm = i-iry-m, Zm = {-ir+' Z_^, m > 0. 



16 



WANG AND YANG 



We derive then the reduced equations as follows: 

dyi r)+ ( 2 t 2 I 2 I 2 I 2\ I io\ 

"^212/1 - 922/1 (a;o + 2/i + + 2/2 + ^2) + o(3) 



dt 

'^^^ -/52I2/2 - q2y2{xl + yl + zl + yi + zl) + 0(3), 



dt 

(5-15) ^ =/3^ia;o - 922:0(2;^ + yi + zl + yl + zl) + o(3), 

dz 

=p}^Zi - q2Zi{xl +yl + zl +yl + z^) + o(3), 

fly 

=PtiZ2 - q2Z2{xl + yl + zl + yl + zl) + o(3). 

Then as in the proof of last theorem, we can show that Sjj is homeomorphic to S^, and we have 
the following approximation: 

Sfl ~ {(?y2,yi,a;o,2;i,2;2) e R'"^ | + yl + zj + yl + zl = [3:^^/q2}, 

and Si{ contains at least an S"^ subset of steady state solutions. We remark here that if we retain 
only the linear and cubic order terms in the above reduced equations (5.15), we would have the whole 
5** consisting only steady state solutions. However, since the steady states solutions are degenerate 
(with zero Jacobian), there is no persistence and stability of these steady state solutions. 
The proof of the theorem is complete. 

□ 

6. Large-Scale Circulations and Turbulent Frictions 
Now we consider the following equations: 



1 du (P' 

^{ul + VuU + w— + Vp) + (Tqu - (A + — 
Pr dz oz' 



(6.1) ^(-. + + + 1) + - + ^> -Vrt = o, 

dw I — 

Tt + V„T + w— - VRw - (A + = 0, 

dw 

divu + —— = 0. 
az 

supplemented with boundary condition: 

(6.2) w = 0,T = 0, — = at z = 0,1. 

oz 

By similar computation, we get the three groups of eigenvalues and eigenvectors: 
(1) For (w,T)=0 

fSio = - Pr{af + CTo), 
*imO = (curiy;„,0,0); 



(2) For / = with n > 1 



*oon = (0,0,sinn7r2) 
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(3) For l,n>l 

PL 



Imn 



Pin 



Here 



= {n-KCOsn-KzWim, OiYimSmmrz, 
--D--\/D'^-AVvE, 
= {n-KCOsn'KzSJYijn, OiYimSmmrz, 

£> = 7?J1 + Pr) + 



2Aaf 



Ptn+lln 



2Xaf 



Pin + -Tin 



Ylm sin 71772;), 



Yi^msmn-Kz). 



iln 



To get the critical Rayleigh number Ac, let /3;,„ = 0, and obtain that 
(6.3) 



A. 



min — \/ (7I + <7iaf + aon^TT^) = min . 



a? 



By the analysis for the convection scales in [13], we consider the case where the turbulent friction 

coefficients satisfy 

(6.4) K CTo < ai. 

In this case, the minimum in (6.3) is approximately achieved at 



(6.5) 



lc{lc + 1) = a' ~ Y 



CTl 



1/2 



Then the same results as Theorems 4.1, 5.1 and 5.2 hold true for the Rayleigh-Bcnard convection 
system (6.1) with added turbulent friction terms. We omit the details. It is worth, however, 
to remark that (6.5) provides a precise pattern selection mechanism in terms of the aspect ratio 
h/a, and the ratio between the horizontal and vertical friction coefficients. For example, for the 
circulation, we can choose the parameters as follows: 

a = 6.4 X lO^m, h = lO^'m, lc = Q 

Then, 

ao/cTi « 10"^ 

Here Ic = 6 represents the 6 Walker circulation cells over the tropical, and is also consistent with 
the three meridional circulation cells-the Hadley cell, the midlatitude cell and the polar cell, in the 
Northern hemisphere. 
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